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ABSTRACT 

Using  imaging  data  from  the  STEREO  mission,  we  empirically  reconstruct 
the  time-dependent  three-dimensional  morphology  of  a  coronal  mass  ejection 
(CME)  from  2008  June  1,  which  exhibits  significant  variation  in  shape  as  it 
travels  from  the  Sun  to  1  AU,  requiring  us  to  abandon  the  assumption  of  self¬ 
similar  expansion.  We  model  the  CME  as  a  flux  rope  that  is  rather  fat  relative  to 
its  longitudinal  extent  close  to  the  Sun,  but  which  becomes  thinner  and  flatter  on 
top  as  the  flux  rope  moves  outwards.  We  find  best  agreement  with  the  STEREO 
images  when  the  flux  rope’s  west  leg  is  assumed  to  be  rotated  35°  below  the 
ecliptic  plane.  This  orientation  is  consistent  with  previously  published  inferences 
about  this  CME’s  orientation  from  an  analysis  of  in  situ  observations  of  the  event 
from  June  6-7,  when  the  CME  hits  STEREO-B.  The  agreement  between  these  two 
very  different  kinds  of  analysis  is  encouraging.  Close  to  1  AU  the  CME  not  only 
hits  STEREO-B  but  also  strikes  a  comet  (Comet  C/2007  W1  Boattini),  which 
provides  an  additional  constraint  for  our  reconstruction  efforts.  Finally,  we  find 
that  this  CME  is  very  instructive  for  assessing  different  methods  of  extracting 
kinematic  information  from  measurements  of  elongation  angles  from  the  Sun, 
which  is  a  complicated  issue  for  measurements  far  from  the  Sun.  The  “Fixed-0” 
assumption  that  we  have  used  successfully  in  the  past  does  not  work  well  here, 
and  we  discuss  the  implications  for  extracting  reliable  kinematic  information  from 
heliospheric  imaging. 

Subject  headings:  Sun:  coronal  mass  ejections  (CMEs)  —  solar  wind  —  inter¬ 
planetary  medium 
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1.  INTRODUCTION 

The  primary  goal  of  the  Solar  TErrestrial  RElations  Observatory  (STEREO)  mission 
(Kaiser  et  al.  2008)  is  to  improve  our  understanding  of  the  three-dimensional  morphology 
of  coronal  mass  ejections  (CMEs).  The  deep,  extended  solar  minimum  of  2008-09  did  not 
provide  many  dramatic  CMEs  for  the  two  STEREO  spacecraft  to  study,  but  on  2008  June  1, 
a  very  slow  coronal  mass  ejection  (CME)  is  observed  drifting  slowly  away  from  the  Sun,  with 
the  event  clearly  discernible  in  images  from  both  spacecraft.  The  CME  expands  off  the  east 
limb  as  viewed  from  STEREO-A,  which  at  the  time  was  located  25°  ahead  of  the  Earth  in 
its  orbit  around  the  Sun.  In  contrast,  the  CME  appears  as  a  halo  event  from  the  perspective 
of  STEREO-B,  orbiting  29°  behind  the  Earth. 

The  imagers  that  make  up  the  SECCHI  instrument  packages  on  each  STEREO  space¬ 
craft  (Howard  et  al.  2008)  successfully  track  this  CME  continously  from  close  to  the  Sun  to 
1  AU.  The  CME  hits  STEREO-B  on  2008  June  6,  where  the  onboard  PLASTIC  (Galvin 
et  al.  2008)  and  IMPACT  (Acuna  et  al.  2008;  Luhmann  et  al.  2008)  instrument  suites  ob¬ 
serve  event  characteristics  typical  of  a  magnetic  flux  rope  (Mostl  et  al.  2009).  The  STEREO 
spacecraft  therefore  furnish  simultaneous  imaging  and  in  situ  observations  of  this  particular 
CME,  providing  a  perfect  opportunity  to  compare  the  appearance  of  a  flux  rope  CME  in 
white  light  images  with  its  manifestation  in  the  in  situ  data. 

One  final  factor  that  has  made  this  event  worthy  of  study  is  the  lack  of  any  discernablc 
surface  activity  associated  with  it.  Both  the  in  situ  detection  at  STEREO-B  and  the  eastward 
direction  of  the  CME  in  STEREO-A  images  demonstrate  that  the  halo  appearance  of  the 
CME  seen  from  STEREO-B’s  perspective  is  due  to  a  front-side  event  and  not  a  back-side 
event.  Nevertheless,  the  EUVI  imagers  on  both  STEREO  spacecraft,  which  view  the  solar 
corona  in  various  EUV  bandpasses,  detect  no  flaring  or  filament  eruption  at  the  time  of 
the  CME’s  inception  (Robbrecht  et  al.  2009),  demonstrating  more  clearly  than  ever  before 
that  CME  initiation  does  not  necessarily  require  any  obvious  surface  driver,  at  least  for  the 
very  slow  “streamer  blowout”  CME  class  to  which  this  event  clearly  belongs  (Howard  et  al. 
2008).  The  purpose  of  this  paper  is  to  study  the  morphological  and  kinematic  evolution  of 
this  CME,  from  its  inception  to  its  arrival  at  STEREO-B,  based  on  an  empirical  analysis  of 
STEREO  images  of  the  event  from  both  the  A  and  B  spacecraft. 

2.  MORPHOLOGY  OF  THE  JUNE  1  CME:  QUALITATIVE  ASSESSMENT 

Each  STEREO  spacecraft  carries  four  white-light  telescopes  that  can  observe  CMEs  at 
different  distances  from  the  Sun  (Howard  et  al.  2008).  There  are  two  coronagraphs,  COR1 
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and  C0R2,  that  observe  the  Sun’s  white  light  corona  at  angular  distances  from  Sun-center 
of  0.37°  —  1.07°  and  0.7°  —  4.2°,  respectively,  corresponding  to  distances  in  the  plane  of  the 
sky  of  1.4  —  4.0  R0  for  COR1  and  2.5  —  15.6  R0  for  COR2.  And  there  are  two  heliospheric 
imagers,  HI1  and  HI2,  that  monitor  the  interplanetary  medium  (IPM)  in  between  the  Sun 
and  Earth,  where  HI1  observes  elongation  angles  from  Sun-center  of  3.9°  —  24.1°  and  HI2 
observes  from  19°  —  89°. 

Figure  1  shows  four  STEREO  images  of  the  2008  June  1  event.  This  CME  is  a  halo 
event  from  STEREO-B’s  perspective.  While  apparent  in  COR2-B  movies,  the  halo  is  quite 
faint  in  the  static  COR2-B  image  in  Figure  1,  so  an  outline  is  provided  to  guide  the  eye.  In 
COR2-A  the  CME  has  a  compact  appearance  with  a  vaguely  circular  top.  Such  a  semblance 
is  commonly  seen  for  CMEs  and  is  often  interpreted  as  representing  the  apex  of  a  flux  rope 
viewed  edge-on  (Chen  et  al.  1997;  Wood  et  al.  1999;  Low  2001;  Cremades  &  Bothmer  2004; 
Krall  2007).  However,  the  HI1-A  appearance  of  the  June  1  CME  is  very  different,  looking 
more  like  a  broad  flux  rope  viewed  closer  to  face-on.  This  visage  is  more  consistent  with 
the  flux  rope  orientation  that  Mostl  et  ah  (2009)  infer  for  the  CME  based  on  the  in  situ 
observations  from  STEREO-B.  They  estimate  a  —45°  tilt  of  the  flux  rope  relative  to  the 
ecliptic  plane,  and  argue  that  this  is  consistent  with  the  HI1-A  and  HI2-A  appearance  of  the 
CME.  The  suggestion  is  made  that  perhaps  the  change  in  appearance  between  COR2-A  and 
HI1-A  is  due  to  a  ~  45°  clockwise  rotation  of  the  flux  rope. 

Although  a  rotation  may  be  consistent  with  STEREO-A’s  view  of  the  event,  it  is  ques¬ 
tionable  whether  it  is  the  best  interpretation  of  STEREO-B’s  images.  To  us,  the  shape  of 
the  COR2-B  halo  is  already  suggestive  of  a  roughly  —45°  tilt  of  the  CME  at  early  times. 
Furthermore,  we  see  no  evidence  for  any  clockwise  rotational  motion  of  the  halo  in  COR1-B, 
COR2-B,  or  HI1-B  movies.  Thus,  our  preferred  interpretation  of  the  event  is  that  the  CME 
flux  rope  possesses  a  tilt  right  from  the  start,  but  this  orientation  is  not  evident  at  early 
times  due  to  the  flux  rope  being  initially  very  fat,  making  the  tilt  angle  difficult  to  infer 
from  STEREO-A’s  perspective,  hr  this  paradigm,  which  will  be  presented  in  detail  in  §4, 
the  change  in  appearance  between  COR2-A  and  HI1-A  is  due  to  a  narrowing  of  the  flux  rope 
and  a  flattening  of  its  top  with  time,  as  the  CME  moves  outwards  from  the  Sun. 

Finally,  there  is  a  comet  visible  in  the  HI2-B  image  in  Figure  1,  which  is  struck  by 
the  CME.  This  is  Comet  C/2007  W1  Boattini.  STEREO  has  previously  observed  other 
comets  interacting  with  the  solar  wind  and  transients  within  it  (Vourlidas  et  al.  2007;  Jia 
et  al.  2009).  There  is  little  evidence  of  a  tail  for  Comet  Boattini  during  the  day  leading 
up  to  the  CME’s  passage  over  the  comet,  at  least  not  in  the  HI2-B  images,  but  a  tail 
appears  and  brightens  starting  at  about  4:10  UT  on  June  6.  The  tail  then  experiences 
a  clear  disconnection  event  at  about  8:10  UT,  and  remains  active  and  variable  for  about 
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2  days  thereafter  before  fading  again.  If  we  associate  the  disconnection  event  at  June  6, 
8:10  UT  with  the  timing  of  the  CME’s  primary  density  pulse  hitting  the  comet,  we  can 
use  our  knowledge  of  the  comet’s  orbit  to  pinpoint  the  CME’s  location  at  8:10  UT,  at  least 
at  one  point  along  the  CME  front.  In  this  way  the  comet  provides  a  useful  constraint  for 
our  morphological  modeling.  In  Heliocentric  Earth  Ecliptic  (HEE)  coordinates,  the  comet’s 
position  at  the  time  of  CME  impact  is  (0.898  AU,  —0.096  AU,  —0.157  AU),  based  on  an 
ephemeris  from  http://ssd.jpl. nasa.gov/sbdb. cgi?sstr=C/2007  Wl. 


3.  KINEMATICS  OF  THE  JUNE  1  CME 

Measuring  CME  velocities  and  trajectories  is  crucial  for  understanding  how  CMEs  prop¬ 
agate  through  the  IPM.  It  is  also  important  for  space  weather  forecasting  purposes  in  order  to 
accurately  predict  if  and  when  a  CME  will  hit  Earth.  For  our  purposes,  we  need  a  kinematic 
model  of  the  June  1  CME  to  properly  scale  our  time-dependent  three-dimensional  (3-D) 
morphological  models.  Both  kinematic  and  morphological  understanding  of  the  event  are 
necessary  to  generate  synthetic  images  that  can  be  compared  directly  with  the  real  images, 
providing  the  ultimate  test  of  both  the  kinematic  and  morphological  aspects  of  the  analysis. 

It  is  much  easier  to  measure  a  CME’s  radial  motion  from  the  Sun  when  the  CME  is 
observed  from  the  side  than  if  the  CME  is  moving  directly  towards  or  away  from  an  observer. 
This  means  that  for  the  June  1  CME  STEREO- A  is  in  a  much  more  useful  location  for 
a  kinematic  analysis  than  STEREO-B.  We  measure  the  elongation  angle  of  the  CME  as  a 
function  of  time  from  all  STEREO- A  images  of  the  event,  from  COR1-A  all  the  way  through 
HI2-A.  In  the  COR2-A  and  HI1-A  panels  of  Figure  1,  plus  signs  indicate  how  we  have  marked 
the  leading  edge  of  the  CME. 

A  kinematic  analysis  requires  computing  real  distances,  r,  from  the  measured  elongation 
angles,  e.  There  are  two  simple  methods  that  have  been  used  in  the  past,  the  “Point-P”  and 
“Fixed-0”  methods  (Kahler  &  Webb  2007;  Howard  et  al.  2007,  2008;  Sheelcy  et  al.  2008; 
Wood  et  ah  2009a).  The  former  assumes  an  extremely  broad  spherical  front  centered  on  the 
Sun,  with 

r  —  d  sine,  (1) 

where  d  is  the  distance  from  the  observer  to  the  Sun,  the  observer  being  STEREO- A  in  this 
case.  The  latter  assumes  a  very  narrow  CME,  with 

d  sin  e 

r  =  •  7  — k,  (2) 

sm  (e  +  <p) 

where  </>  is  the  angle  between  the  CME  trajectory  and  the  observer’s  line  of  sight  to  the 
Sun.  More  recently,  Lugaz  et  ah  (2009)  have  proposed  a  third  method,  which  approximates 
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the  CME  as  a  sphere  centered  halfway  between  the  Sun  and  the  CME’s  leading  edge.  This 
assumption  yields 

2d  sine 

r  =  1  +  sin(e  +  d)  ‘  (3) 

Lugaz  et  al.  (2009)  call  this  the  “Harmonic  Mean”  approximation,  since  an  r  value  computed 
under  this  assumption  represents  the  intermediate  harmonic  mean  of  values  computed  under 
the  other  two  approximations. 

The  top  panel  of  Figure  2  shows  the  leading  edge  distances  computed  for  the  June  1 
CME  for  all  three  approximations,  and  the  bottom  panel  shows  the  inferred  velocity  of  the 
CME’s  leading  edge.  The  morphological  analysis  presented  below  in  §4  suggests  a  CME 
trajectory  of  <f>  =  63°,  so  that  is  what  is  assumed  for  the  Fixed-d  and  Harmonic  Mean 
assumptions.  Velocities  are  generally  computed  from  adjacent  distance  data  points,  but 
this  sometimes  leads  to  velocities  with  huge  error  bars,  which  vary  wildly  in  time  in  a  very 
misleading  fashion.  For  this  reason,  as  we  compute  the  velocities  from  left  to  right  in  the 
figure  we  actually  skip  distance  points  until  the  uncertainty  in  the  computed  velocity  ends  up 
under  some  assumed  threshold  value  (100  km  s-1  in  this  case),  similar  to  what  we  have  done 
in  the  past  (Wood  et  al.  1999,  2009a).  The  velocity  uncertainties  are  computed  assuming 
the  following  estimates  for  the  uncertainties  in  the  distance  measurements:  1%  fractional 
errors  for  the  COR1  and  COR2  distances,  and  2%  and  3%  uncertainties  for  HI1  and  H12, 
respectively. 

The  three  velocity  profiles  in  Figure  2  are  consistent  throughout  the  COR1,  COR2,  and 
HI1  fields  of  view,  suggesting  a  CME  that  drifts  through  COR1-A  at  a  roughly  constant 
speed  of  about  50  km  s'1  and  then  accelerates  through  COR2-A  into  HI1-A,  reaching  a 
velocity  of  about  400  km  s-1.  However,  the  three  assumptions  represented  by  equations 
(l)-(3)  result  in  wildly  different  conclusions  about  what  happens  to  the  CME  in  the  HI2-A 
held  of  view.  The  Point-P  approximation  implies  a  dramatic  deceleration  in  H12-A,  while  the 
Fixed-d  approximation  implies  a  dramatic  acceleration.  It  is  very  difficult  to  imagine  any 
physical  reason  why  this  CME  would  do  either  this  far  from  the  Sun.  The  Harmonic  Mean 
approximation  yields  a  far  more  plausible  kinematic  profile,  with  the  CME  speed  simply 
leveling  off  nicely  at  400  km  s-1. 

Aside  from  plausibility  arguments,  another  way  to  validate  the  Harmonic  Mean  approx¬ 
imation’s  distance  and  velocity  measurements  is  to  compare  them  with  the  in  situ  measure¬ 
ments  of  the  CME  by  STEREO-B  (see  Mostl  et  al.  2009).  We  find  that  the  CME  arrival 
time  at  STEREO-B  suggested  by  the  Harmonic  Mean  measurements  is  within  an  hour  of  the 
observed  June  6,  UT  21:00  arrival  time  of  the  CME’s  main  density  peak.  The  1  AU  velocity 
suggested  by  the  Harmonic  Mean  approximation  is  also  consistent  with  the  observed  ~  400 
km  s'1  velocity  at  STEREO-B.  In  contrast,  the  Fixe  d-d  CME  kinematics  have  the  CME 


-6  - 

arriving  about  a  half  day  too  early,  and  at  too  high  a  speed.  The  Point-P  kinematics  have 
the  CME  arriving  way  too  late,  and  at  too  low  a  speed. 

We  have  analyzed  the  kinematics  of  several  CMEs  in  the  past  for  which  the  Fixed-0 
approximation  did  not  lead  to  implausible  velocity  profiles,  and  which  accurately  predicted 
1  AU  arrival  times  (Wood  et  al.  2009a, b;  Wood  &  Howard  2009).  So  what  is  different  about 
this  CME?  Is  it  possible  that  we  are  assuming  the  wrong  CME  trajectory,  0?  After  all,  one 
disadvantage  of  both  the  Fixed-0  and  Harmonic  Mean  methods  is  that  they  require  that  q t> 
be  precisely  known,  and  be  applicable  at  all  times  for  a  given  CME. 

Figure  3  shows  distance  and  velocity  versus  time  computations  for  the  June  1  CME 
assuming  two  different  </>  values,  under  the  Fixed- 0  assumption.  We  already  noted  in  Figure  2 
that  the  <fi  =  63°  Fixed-0  model  yields  an  improbable  HI2-A  acceleration,  and  that  this  can 
be  solved  by  switching  to  the  Harmonic  Mean  approach.  Figure  3  shows  that  another  solution 
is  to  reduce  (j)  by  only  13°  to  </>  =  50°.  One  approach  to  any  kinematic  analysis  is  to  actually 
assume  that  the  correct  trajectory  is  whatever  (j)  leads  to  a  constant  velocity  far  from  the 
Sun,  in  effect  assuming  that  the  CME  has  reached  some  sort  of  a  terminal  velocity  by  the 
time  it  reaches  the  HI2  field  of  view  (Sheeley  et  al.  2008;  Rouillard  et  al  2008).  However,  we 
would  argue  that  for  this  CME  the  trajectory  estimates  close  to  the  Sun  are  quite  robust. 
Our  preferred  (j)  =  63°  comes  from  our  own  morphological  analysis  in  §4,  but  this  result  is 
quite  consistent  with  the  (j)  =  66°  value  that  Thernisien  et  al.  (2009)  infer  from  COR2  data 
alone.  It  is  also  consistent  with  other  assessments  close  to  the  Sun,  which  are  summarized 
by  Mostl  et  al.  (2009). 

Is  it  possible  that  the  CME  trajectory  actually  changes  from  0  63°  close  to  the  Sun 

to  (j)  ~  50°  in  the  HI2  field  of  view?  In  some  sense,  we  believe  that  this  is  indeed  what  is 
happening.  However,  it  is  not  that  the  CME  is  changing  direction,  it  is  instead  that  the 
part  of  the  CME  front  seen  as  the  leading  edge  is  changing  as  the  CME  moves  outwards.  A 
fundamental  problem  with  the  Fixed-0  approach  is  that  it  assumes  that  a  fixed  observer  sees 
the  same  part  of  the  CME  front  as  the  leading  edge  at  all  times.  The  two  panels  of  Figure  4 
illustrate  how  this  can  be  a  problem.  In  the  first  panel  it  is  the  far  side  of  the  CME  front 
that  is  seen  as  the  leading  edge.  This  apparent  leading  edge  shifts  to  the  foreground  side  of 
the  CME  front  as  the  CME  moves  outwards.  We  believe  this  is  why  the  HI2-A  data  seem 
to  suggest  a  lower  trajectory  angle  than  our  preferred  (j)  =  63°  angle.  The  part  of  the  CME 
front  that  is  being  followed  in  HI2-A  truly  is  at  a  lower  trajectory  angle  than  the  actual 
center  of  the  CME. 

The  Point-P  and  Harmonic  Mean  approximations  both  allow  the  apparent  leading  edge 
of  the  CME  to  shift  along  the  CME  front,  albeit  with  specific  assumptions  about  the  shape 
of  the  front.  For  this  particular  CME,  the  allowances  provided  by  the  Harmonic  Mean 
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assumption  seem  to  have  led  to  a  proper  reproduction  of  the  actual  CME  kinematic  behavior. 
However,  none  of  the  three  techniques  represented  by  equations  (l)-(3)  should  be  considered 
superior.  Each  technique  will  be  applicable  for  particular  viewing  geometries  and  certain 
CME  shapes  (e.g.,  Howard  &  Tappin  2009).  The  Fixed-0  approach  will  be  best  for  narrow 
CMEs,  or  cases  in  which  only  a  small  part  of  a  broad  CME  front  is  visible.  The  Point-P 
and  Harmonic  Mean  approximations  will  be  more  applicable  for  broad,  bright  CME  fronts, 
with  the  Point-P  approximation  probably  being  better  for  CMEs  directed  more  towards  the 
observer.  A  plot  like  Figure  2  should  be  helpful  in  assessing  which  technique  is  best  for  a 
given  event,  especially  if  the  assumption  is  made  that  the  best  kinematic  model  is  that  which 
yields  the  flattest  velocity  profile  far  from  the  Sun.  If  the  CME  hits  a  spacecraft,  the  arrival 
time  information  provided  by  in  situ  data  will  also  be  helpful. 

In  Figure  5,  we  have  fitted  a  kinematic  model  to  the  Harmonic  Mean  distance-vs.-time 
data  points  from  Figure  2.  This  is  a  three-phase  model,  somewhat  analogous  to  kinematic 
models  we  have  used  to  model  CMEs  in  the  past  (Wood  et  al.  2009a, b),  but  with  differences 
associated  with  the  very  slow  initial  speed  of  the  CME.  The  model  assumes  an  initial  phase 
of  constant  velocity,  V\,  ending  at  some  time  7 j ;  a  second  phase  of  constant  acceleration, 
a2,  ending  at  some  time  t-2]  followed  by  a  final  phase  of  constant  velocity.  In  addition  to 
these  four  free  parameters,  there  is  a  fifth  involving  a  time  shift  from  the  model  time  vector 
to  that  of  the  actual  measurements.  Figure  5  shows  the  best  fit  to  the  data  based  on  this 
simple  five-parameter  model.  The  CME  accelerates  from  35  km  s-1  at  a  Sun-center  distance 
of  3.7  R0  to  384  km  s-1  at  28  R0,  and  then  maintains  that  speed  all  the  way  out  to  1  AU, 
consistent  with  the  point-by-point  kinematic  behavior  seen  in  Figure  2. 


4.  MORPHOLOGY  OF  THE  JUNE  1  CME:  A  3-D  RECONSTRUCTION 


Our  primary  goal  here  is  to  empirically  derive  a  model  of  the  3-D  mass  distribution 
of  the  CME  assuming  a  flux  rope  topology,  a  model  which  is  designed  to  reproduce  the 
CME’s  appearance  in  the  STEREO  images  as  well  as  possible.  For  this  purpose,  we  use 
a  mathematical  formulation  of  a  flux  rope  shape  that  we  have  used  in  the  past  (Wood  et 
al.  2009a).  The  shape  of  the  inner  and  outer  edges  of  a  2-D  flux  rope  are  defined  in  polar 
coordinates  using  the  following  equation: 


r(9) 


=  rmax  exp 


(4) 


These  two  2-D  loops  are  used  to  define  a  3-D  flux  rope  shape  by  assuming  a  circular  cross 
section  for  the  flux  rope,  bounded  by  the  two  loops. 


The  flux  rope  shape  is  then  converted  to  Cartesian  coordinates  and  then  densities  are 


mapped  into  a  density  cube  by  assuming  a  Gaussian  density  profile  across  the  flux  rope 
surface: 


ni(x,y,z)  =  nmax  exp 


5(x,y,z) 


(Jr 


(5) 


where  <5(x,  y ,  z)  is  the  distance  to  the  flux  rope  surface  from  a  point  (x,y,z)  within  the  density 
cube.  To  allow  for  an  x-dependence  for  the  density,  the  density  cube  is  further  modified  by 


n2(x,y,z)  =  m(x,y,z) 


x 


(B) 


The  final  step  is  to  rotate  the  flux  rope  to  the  desired  direction.  This  means  first  defining  a 
longitude  and  latitude  (lF,  bF)  that  indicate  the  direction  of  the  CME,  or  more  specifically 
the  direction  of  the  x-axis  of  the  density  cube.  Because  the  flux  rope  is  not  rotationally 
symmetric,  the  flux  rope  requires  an  additional  angle,  <j)fr ,  defining  the  orientation  of  the 
flux  rope  with  respect  to  the  ecliptic  plane.  Once  a  density  cube  with  a  3-D  model  flux  rope 
is  complete,  we  generate  synthetic  images  from  it  using  a  white-light  rendering  routine  to 
perform  the  necessary  calculations  of  Thomson  scattering  within  the  density  cube  (Billings 
1966;  Thcrnisien  et  al.  2006;  Wood  et  al.  2009b).  The  various  parameters  of  the  flux  rope 
are  adjusted  to  maximize  the  agreement  between  the  real  and  synthetic  images  of  the  CME. 


The  above  functional  form  for  flux  rope  CMEs  was  first  used  to  model  the  flux  rope 
shape  in  a  CME  from  2008  April  26  (Wood  et  ah  2009a).  In  that  case,  we  found  that  the 
appearance  of  the  CME  could  be  reproduced  reasonably  well  by  assuming  that  the  flux  rope 
simply  expands  in  a  self-similar  fashion.  In  reality,  there  was  evidence  for  some  evolution 
in  the  CME’s  morphology  with  time,  but  the  self-similar  expansion  assumption  was  at  least 
a  decent  first-order  approximation.  Based  on  Figure  1  it  is  immediately  apparent  that  this 
is  not  the  case  for  the  2008  June  1  CME.  For  this  CME,  the  CME  shape  must  be  allowed 
to  evolve  with  time  by  allowing  some  of  the  model  parameters  to  vary,  complicating  the 
analysis. 


Table  1  lists  the  model  parameters  that  lead  to  our  best  reproduction  of  the  CME’s 
appearance  in  the  STEREO  images,  as  judged  by  eye.  The  (lF,bF)—( 217°, 2°)  angles  indicate 
the  3-D  CME  trajectory  in  heliocentric  aries  ecliptic  (HAE)  coordinates.  This  corresponds  to 
a  direction  63°  from  STEREO-A,  indicating  where  the  =  63°  assumption  in  §3  originates. 
It  is  the  kinematic  model  of  the  CME  described  in  §3  that  establishes  the  absolute  size  scale 
of  our  density  cube  as  a  function  of  time  (specifically  the  kinematic  model  in  Fig.  5),  so  the 
rma x  of  the  outer  edge  is  simply  normalized  to  1  and  other  distances  in  the  parameter  list 
are  quoted  relative  to  this  value.  Since  we  are  only  trying  to  reproduce  the  CME’s  observed 
morphology  and  kinematic  behavior,  and  not  its  absolute  brightness,  the  absolute  value  of 
the  density  is  not  of  interest,  so  nmax  is  simply  set  to  1. 
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To  provide  for  the  required  evolution  of  the  flux  rope  shape,  there  are  three  parameters 
in  Table  1  that  vary  with  time,  and  ranges  ([aq,^])  are  quoted  for  these  parameters  in  the 
table.  The  rmax  value  of  the  flux  rope’s  inner  edge  increases  with  time,  resulting  in  a  flux 
rope  that  gradually  becomes  thinner  relative  to  its  length.  The  a  of  the  outer  edge  increases 
with  time,  yielding  a  flux  rope  that  flattens  at  its  top.  Finally,  the  a  of  the  outer  edge  is 
allowed  to  increase  slightly,  meaning  the  longitudinal  extent  of  the  flux  rope  increases  slightly 
with  time. 

The  following  equation  describes  how  the  time-dependent  parameters  are  computed 
from  the  [aq,2q]  ranges  quoted  in  the  table: 

x(t)  =  aq  +  X 2  -  Xl  tanh  +  1  •  (7) 

The  reference  time,  t0,  and  relaxation  time,  tr,  are  additional  fit  parameters,  which  are  also 
listed  in  Table  1.  Assuming  t  —  0  corresponds  to  11:45  UT  on  2008  June  1,  we  settle  on 
values  of  to  =  40  hr  and  tr  =  15  hr. 

Figure  6  compares  real  and  synthetic  images  at  three  separate  times  in  the  COR2, 
HI1,  and  HI2  fields  of  view,  respectively.  Both  the  real  and  synthetic  images  are  displayed 
in  running  difference  format,  where  the  previous  image  is  subtracted  from  each  image  to 
emphasize  the  dynamic  CME  structure  as  much  as  possible.  Individual  panels  also  illustrate 
the  3-D  shape  of  the  model  structure,  the  position  of  the  spacecraft  and  model  CME  in  the 
ecliptic  plane,  and  observed  solar  wind  properties  at  STEREO-A  and  STEREO-B  at  the 
time  of  the  images.  The  data-modcl  comparison  is  fairly  straightforward  for  the  STEREO-A 
images,  but  the  faintness  of  the  CME  in  STEREO-B  images  makes  it  harder  to  assess  the 
quality  of  fit  to  those  data.  For  STEREO-B  in  particular,  movies  are  greatly  preferable  to 
static  images  in  comparing  the  real  and  synthetic  data,  so  a  movie  version  of  Figure  6  has 
been  provided  in  the  online  version  of  this  article. 

The  COR2  frame  of  Figure  6  illustrates  how  the  CME’s  appearance  at  early  times  can 
be  reproduced  by  a  flux  rope  with  a  significant  35°  tilt  relative  to  the  ecliptic,  as  long  as 
it  is  a  fat  flux  rope.  The  two  subsequent  HI1  and  HI2  frames  illustrate  how  narrowing 
and  flattening  the  flux  rope  can  explain  the  evolving  shape  of  the  CME  in  the  STEREO-A 
images,  although  only  the  central  part  of  the  CME  is  really  visible  in  the  HI2-A  image.  This 
demonstrates  that  it  is  possible  to  reproduce  the  changing  appearance  of  the  CME  without 
requiring  any  rotation  of  the  structure  between  the  COR2-A  and  HI1-A  fields  of  view.  The 
HI2  frame  also  shows  the  location  of  the  3-D  flux  rope  relative  to  Earth,  the  two  STEREO 
spacecraft,  and  Comet  Boattini  at  a  time  shortly  before  the  CME  strikes  STEREO-B  and 
the  comet. 

Although  the  model  reproduces  the  shape  of  the  outer  edge  of  the  CME  flux  rope 
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reasonably  well,  the  HI1  panel  of  Figure  6  shows  that  the  model  is  not  as  successful  in  fitting 
the  shape  and  latitudinal  extent  of  the  inner  edge.  The  outline  of  the  flux  rope  inner  edge 
suggested  by  the  model  is  much  too  pointed.  This  is  an  indication  of  the  limitations  of  our 
parametrized  flux  rope  shape,  as  no  combination  of  parameters  will  rectify  this  problem.  To 
be  specific,  equation  (4)  simply  will  not  yield  an  inner  edge  shape  that  will  better  reproduce 
the  observations.  A  more  precise  replication  would  probably  require  a  functional  form  that 
allows  the  legs  of  the  loop  to  start  very  close  together  and  then  expand  away  from  each  other 
superradially,  something  equation  (4)  does  not  allow. 

The  STEREO  mission  promises  to  greatly  improve  our  understanding  of  the  angular 
extent  of  CMEs  in  the  inner  heliosphere.  An  angular  extent  for  our  model  CME  in  the  plane 
of  the  flux  rope  can  be  estimated  by  first  setting  r(6)  =  0.9 rmax  in  equation  (4),  and  then 
solving  for  9,  assuming  the  outer  edge  o  and  a  parameters  listed  in  Table  1.  The  resulting 
9  represents  the  half-width  of  the  flux  rope,  so  it  should  be  multiplied  by  2  to  yield  the  full 
width.  This  exercise  leads  to  an  angular  width  estimate  of  70°  for  the  June  1  CME.  This  can 
be  compared  with  the  angular  extents  inferred  for  various  previous  events  based  on  in  situ 
observations.  Bothmer  &  Schwenn  (1998)  find  that  magnetic  clouds  observed  by  the  Helios 
spacecraft  are  generally  also  observed  at  Earth  if  less  than  60°  away.  Richardson  &  Cane 
(1993)  infer  angular  extents  of  ~  100°  for  many  interplanetary  shocks.  Our  70°  flux  rope 
might  at  first  seem  consistent  with  these  numbers,  but  it  is  important  to  note  that  due  to 
the  tilt  of  the  June  1  CME  flux  rope,  in  the  ecliptic  plane  the  CME  only  has  an  extent  of 
~  25°  (see  panel  F  of  the  HI2  frame  in  Fig.  6),  explaining  why  the  event  is  not  seen  at  Earth 
only  29°  from  STEREO-B.  Thus,  we  suspect  the  generally  more  energetic  events  studied  by 
Bothmer  &  Schwenn  (1998)  and  Richardson  &  Cane  (1993)  are  probably  larger  on  average 
than  the  2008  June  1  CME. 

The  in  situ  data  panel  in  the  H12  frame  shows  the  two  density  peaks  seen  by  STEREO- 
B  when  the  CME  hits,  which  are  discussed  in  much  more  detail  by  Mostl  et  al.  (2009). 
The  two  peaks  correspond  nicely  to  the  two  parallel  fronts  observed  in  the  H12-A  image, 
representing  the  front  and  back  side  of  the  flux  rope.  In  movies,  two  broad  but  faint  fronts 
are  also  observed  to  pass  through  the  HI2-B  held  of  view  corresponding  to  these  two  density 
peaks,  at  the  time  when  STEREO-B  is  being  hit  by  the  CME  (see  HI2-B  panel  of  Fig.  1). 
This  illustrates  the  HI2  telescope’s  unprecedented  ability  to  image  a  flux  rope  CME  from 
the  inside.  The  synthetic  images  from  the  model  are  at  least  crudely  able  to  match  the 
appearance  and  timing  of  the  HI2-B  fronts,  the  first  of  which  is  most  evident  in  its  effect  on 
Comet  Boattini. 


According  to  the  model,  the  comet  is  struck  by  the  west  leg  of  the  flux  rope,  while 
STEREO-B  is  struck  by  the  central  part,  with  the  apex  of  the  flux  rope  passing  a  little 
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to  the  south  of  the  spacecraft.  Our  model  CME  first  hits  STEREO-B  at  about  22:10  UT 
on  June  6,  within  an  hour  of  the  CME  front’s  arrival  based  on  the  in  situ  data.  The 
timing  of  the  CME’s  passage  over  the  comet  is  not  quite  as  successfully  reproduced,  with 
the  model  CME  hitting  the  comet  about  4  hours  after  the  observed  tail  disconnection  event. 
Nevertheless,  requiring  that  our  flux  rope  model  strike  the  comet  at  roughly  the  appropriate 
time  represents  a  useful  constraint  on  the  model.  Images  at  early  times  suggest  to  us  a  flux 
rope  tilt  of  (f)fr  ~  —45°,  but  in  our  best  fit  we  reduce  the  tilt  to  <f>fr  =  —35°  to  keep  the 
CME’s  west  leg  from  being  too  much  to  the  south  to  hit  the  comet  squarely. 

Mostl  et  al.  (2009)  estimate  a  flux  rope  tilt  of  about  4>fr  =  —45°  based  on  an  analysis 
of  STEREO-B ’s  in  situ  measurements  of  the  event.  Our  assessment  solely  from  the  images 
is  nicely  consistent  with  their  measurements.  It  is  very  encouraging  that  two  very  different 
methods  of  analyzing  CME  flux  rope  orientation  should  yield  compatible  results. 

Most  of  the  morphological  evolution  of  this  CME  seems  to  occur  in  the  HI1  held  of 
view.  It  is  probably  not  a  coincidence  that  this  is  also  where  the  CME  reaches  its  terminal 
velocity  of  about  400  km  s-1.  The  morphological  changes  are  presumably  driven  mostly  by 
the  CME’s  interactions  with  the  ambient  solar  wind,  but  perhaps  internal  forces  within  the 
CME’s  magnetic  structure  could  also  play  a  role.  Reconstructions  of  CME  evolution  like 
that  performed  here  could  potentially  provide  modelers  with  valuable  information  in  their 
efforts  to  understand  how  CMEs  expand  into  the  IPM  and  interact  with  their  environment. 

In  §3,  we  discussed  at  length  how  different  assumptions  about  the  shape  of  the  leading 
edge  of  the  CME  can  lead  to  very  different  inferences  about  its  kinematic  behavior.  The 
detailed  morphological  reconstruction  of  the  CME  in  this  section  should  in  principle  allow 
an  even  more  precise  connection  to  be  established  between  measured  elongation  angle,  e, 
and  inferred  Sun-center  distance,  r.  In  other  words,  with  the  results  of  the  morphological 
analysis  in  hand,  we  should  be  able  to  do  better  than  the  “Harmonic  Mean”  approximation 
that  we  decided  in  §3  was  the  best  assumption  for  this  particular  CME. 

Using  experiments  such  as  that  in  Figure  4  (but  this  time  using  the  real  leading  edge 
shape  from  our  morphological  analysis),  we  can  establish  exactly  what  part  of  the  CME  front 
is  being  perceived  as  the  leading  edge  as  a  function  of  time,  and  we  can  therefore  connect 
the  apparent  leading  edge  distance  to  the  actual  distance  of  the  central  part  of  the  CME 
front.  Figure  7  shows  that  correcting  our  leading  edge  distances  in  this  fashion  leads  to  a 
kinematic  behavior  far  from  the  Sun  that  is  certainly  more  plausible  than  the  simple  Fixed-0 
and  Point-P  approaches,  and  is  similar  to  the  preferred  Harmonic  Mean  results  (see  Fig.  2). 
However,  despite  the  sophistication  of  this  approach  to  the  kinematic  analysis,  these  new 
distance-vs.-time  measurements  are  not  clearly  superior  to  those  from  the  Harmonic  Mean 
approximation.  In  fact,  a  more  detailed  kinematic  analysis  to  the  new  distance  measurements 
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in  Figure  7  (analogous  to  that  in  Figure  5)  yields  a  velocity  and  arrival  time  prediction  at 
STEREO-B  that  are  in  slightly  worse  agreement  with  the  STEREO-B  in  situ  observations 
than  the  predictions  from  the  Harmonic  Mean  approximation. 

We  conclude  that  the  complicated  morphological  corrections  to  CME  kinematics  repre¬ 
sented  by  Figure  7  may  not  yield  better  results  than  the  simpler  approximations  shown  in 
Figure  2,  and  therefore  will  not  generally  be  worth  the  effort.  The  morphological  corrections 
carry  their  own  questionable  assumptions,  despite  the  sophistication  of  the  approach.  In 
particular,  the  analysis  implicitly  assumes  that  the  entire  CME  leading  edge  is  visible  at  all 
times.  In  reality,  the  HI2-A  images  of  the  CME  tell  us  that  the  most  foreground  parts  of  the 
CME  fade  as  the  CME  moves  through  the  HI2-A  field  of  view,  in  contrast  to  the  synthetic 
images.  Presumably,  there  is  simply  not  enough  mass  there  for  those  parts  of  the  CME 
front  to  be  visible  as  the  CME  moves  farther  from  the  Sun.  This  introduces  errors  in  our 
inferences  about  what  part  of  the  CME  front  is  being  seen  as  the  leading  edge  as  a  function 
of  time.  And  we  believe  these  errors  are  primarily  responsible  for  the  morphological  analysis 
not  leading  to  better  kinematics  compared  to  the  Harmonic  Mean  approximation. 


5.  SUMMARY 

We  have  analyzed  STEREO  observations  of  the  changing  morphology  of  the  2008  June  1 
CME,  and  come  to  the  following  conclusions: 

1.  After  expanding  at  a  constant  slow  velocity  of  about  35  km  s-1  through  the  COR1-A 
field  of  view,  the  leading  edge  of  the  CME  accelerates  through  COR2-A  to  a  terminal 
velocity  of  about  384  km  s^1  28  R0  from  Sun-center  in  the  HI1-A  held  of  view. 

2.  The  observed  change  in  morphology  is  best  explained  by  a  flux  rope  rotated  —35° 
from  the  ecliptic  that  flattens  and  narrows  with  time,  with  the  greatest  changes  taking 
place  at  about  the  time  the  CME  reaches  its  terminal  velocity.  No  change  in  flux  rope 
orientation  is  required  to  explain  the  data,  only  a  change  in  its  shape.  The  CME  has 
an  angular  extent  of  about  70°  in  the  plane  of  the  flux  rope. 

3.  Our  inferred  flux  rope  orientation  agrees  well  with  that  inferred  by  Mostl  et  al.  (2009) 
based  solely  on  in  situ  observations  of  the  CME  from  STEREO-B. 

4.  We  infer  a  CME  trajectory  of  <f>  =  63°  from  STEREO-A’s  location,  consistent  with 
previously  published  assessments. 

5.  The  morphology  and  kinematics  of  our  reconstructed  CME  reproduces  the  timing  of 
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the  CME’s  observed  interaction  with  Comet  Boattini  to  within  4  hours,  and  the  timing 
of  its  arrival  at  STEREO-B  to  within  an  hour. 

6.  Using  the  Fixed-0  approximation  to  infer  the  kinematics  of  this  CME  leads  to  the 
inference  of  an  implausible  acceleration  far  from  the  Sun.  In  contrast,  the  Harmonic 
Mean  approximation  proposed  by  Lugaz  et  al.  (2009)  leads  to  more  plausible  kine¬ 
matics,  and  successfully  reproduces  the  CME’s  velocity  and  arrival  time  observed  by 
STEREO-B  at  1  AU. 

7.  Although  the  Fixed-0  approximation  does  not  work  well  for  this  event,  it  has  worked 
well  for  other  events  that  we  have  analyzed.  The  implication  is  that  there  is  no  single 
best  method  for  getting  from  measured  elongation  angles  to  inferred  CME  Sun-center 
distances.  Many  methods  should  be  tried  for  each  event,  with  the  method  that  leads 
to  a  constant  velocity  far  from  the  Sun  presumably  being  the  best  approximation. 


The  STEREO/SECCHI  data  are  produced  by  a  consortium  of  NRL  (US),  LMSAL  (US), 
NASA/GSFC  (US),  RAL  (UK),  UBHAM  (UK),  MPS  (Germany),  CSL  (Belgium),  IOTA 
(France),  and  IAS  (France).  In  addition  to  funding  by  NASA,  NRL  also  received  support 
from  the  USAF  Space  Test  Program  and  ONR.  In  addition  to  SECCHI,  this  work  has  made 
use  of  data  provided  by  the  STEREO  PLASTIC  and  IMPACT  teams,  supported  by  NASA 
contracts  NAS5-00132  and  NAS5-00133.  We  have  also  made  use  of  data  provided  by  the 
CELIAS/MTOF  experiment  on  SOHO,  which  is  a  joint  ESA  and  NASA  mission.  We  thank 
the  ACE  SWEPAM  and  MAG  instrument  teams  and  the  ACE  Science  Center  for  providing 
the  ACE  data. 


REFERENCES 

Acuna,  M.  H.,  Curtis,  D.,  Scheifele,  J.  L.,  Russell,  C.  T.,  Schroeder,  P.,  Szabo,  A.,  & 
Luhmann,  J.  G.  2008,  Space  Sci.  Rev.,  136,  203 

Billings,  D.  E.  1966.  A  Guide  to  the  Solar  Corona  (New  York:  Academic  Press) 

Bothmer,  V.,  &  Schwenn,  R.  1998,  Ann.  Geophys.,  16,  1 

Chen,  J.,  et  al.  1997,  ApJ,  490,  L191 

Cremades,  H.,  &  Bothmer,  V.  2004,  A&A,  422,  307 

Galvin,  A.  B.,  et  al.  2008,  Space  Sci.  Rev.,  136,  437 


-14 

Howard,  R.  A,  et  al.  2008,  Space  Sci.  Rev.,  136,  67 

Howard,  R.  A.,  Sheeley,  N.  R.,  Jr.,  Michels,  D.  J.,  &  Koomen,  M.  J.  1985,  J.  Geophys.  Res., 
90,  8173 

Howard,  T.  A.,  Fry,  C.  D.,  Johnston,  J.  C.,  &  Webb,  D.  F.  2007,  ApJ,  667,  610 
Howard,  T.  A.,  Nandy,  D.,  &  Koepke,  A.  C.  2008,  J.  Geophys.  Res.,  113,  A01104 
Howard,  T.  A.,  &  Tappin,  S.  J.  2009,  Space  Sci.  Rev.,  147,  31 
Jia,  Y.  D.,  et  al.  2009,  ApJ,  696,  L56 

Kahler,  S.  W.,  &  Webb,  D.  F.  2007,  J.  Geophys.  Res.,  112,  A09103 

Kaiser,  M.  L.,  Kucera,  T.  A.,  Davila,  J.  M.,  St.  Cyr,  O.  C.,  Guhathakurta,  M.,  &  Christian, 
E.  2008,  Space  Sci.  Rev.,  136,  5 

Krall,  J.  2007,  ApJ,  657,  559 

Low,  B.  C.  2001,  J.  Geophys.  Res.,  106,  25141 

Lugaz,  N.,  Vourlidas,  A.,  &  Roussev,  I.  I.  2009,  Ann.  Geophys.,  27,  3479 
Luhmann,  J.  G.,  et  al.  2008,  Space  Sci.  Rev.,  136,  117 

Mostl,  C.,  Farrugia,  C.  J.,  Temmer,  M.,  Miklenic,  C.,  Veronig,  A.  M.,  Galvin,  A.,  Leitner, 
M.,  &  Biernat,  H.  K.  2009,  ApJ,  705,  L180 

Richardson,  I.  G.,  &  Cane,  H.  V.  1993,  J.  Geophys.  Res.,  98,  15295 

Robbrecht,  E.,  Patsourakos,  S.,  &  Vourlidas,  A.  2009,  ApJ,  701,  283 

Rouillard,  A.  P.,  et  al.  2008,  Geophys.  Res.  Lett.,  35,  L10110 

Sheeley,  N.  R.,  Jr.,  et  al.  2008,  ApJ,  675,  853 

Thernisien,  A.  F.  R.,  Howard,  R.  A.,  &  Vourlidas,  A.  2006,  ApJ,  652,  763 

Thernisien,  A.,  Vourlidas,  A.,  &  Howard,  R.  A.  2009,  Sol.  Phys.,  256,  111 

Vourlidas,  A.,  Davis,  C.  J.,  Eylcs,  C.  J.,  Crothers,  S.  R.,  Harrison,  R.  A.,  Howard,  R.  A., 
Moses,  J.  D.,  &  Socker,  D.  G.  2007,  ApJ,  668,  L79 

Wood,  B.  E.,  &  Howard,  R.  A.  2009,  ApJ,  702,  901 


15 


Wood,  B.  E.,  Howard,  R.  A.,  Plunkett,  S.  P.,  &  Socker,  D.  G.  2009a,  ApJ,  694,  707 

Wood,  B.  E.,  Howard,  R.  A.,  Thernisien,  A.,  Plunkett,  S.  P.,  &  Socker,  D.  G.  2009b,  Sol. 
Phys.,  259,  163 

Wood,  B.  E.,  Karovska,  M.,  Chen,  J.,  Brneckner,  G.  E.,  Cook,  J.  W.,  &  Howard,  R.  A.  1999, 
ApJ,  512,  484 


This  preprint  was  prepared  with  the  AAS  PTgX  macros  v5.2. 


16 


Table  1.  CME  Morphological  Parameters 


Quantity 

Flux  Rope 

(outer  edge)  (inner  edge) 

T max 

1  [0.55,0.84]a 

a  (deg) 

[38.2,42.5]a  21.2 

a 

[2.0,8.0]a  4.0 

Tl max 

1 

0.034 

P 

3 

Pfr  (deg) 

-35 

lF  (deg) 

217 

bF  (deg) 

2 

to  (hr) 

40 

tr  (hr) 

15 

aRanges  are  provided  for  time- 
dependent  parameters  (see  text). 
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Fig.  1. —  Four  sample  images  of  the  2008  June  1  CME,  the  left  two  from  STEREO-A  and 
the  right  two  from  STEREO-B.  Note  the  dramatic  change  in  CME  appearance  between  the 
COR2-A  and  HI1-A  images.  The  CME  is  a  faint  halo  event  from  the  perspective  of  STEREO- 
B,  and  a  dotted  line  is  used  to  mark  the  faint  halo’s  location  in  the  COR2-B  image.  The 
HI2-B  image  shows  the  broad,  and  very  faint  CME  front  hitting  Comet  Boattini,  affecting 
the  behavior  of  the  comet’s  tail.  The  positions  of  Earth  and  STEREO-A  are  also  indicated 
in  the  H12-B  image. 
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Fig.  2. —  Distance  and  velocity  measurements  for  the  2008  June  1  CME’s  leading  edge, 
for  three  different  approaches  to  extracting  distances  from  elongation  angle  measurements. 
The  “Fixed-0”  method  yields  an  implausible  acceleration  near  1  AU,  while  the  “Point-P” 
method  produces  an  unlikely  deceleration.  The  kinematics  suggested  by  the  “Harmonic 
Mean”  approach  are  more  plausible. 
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Fig.  3. —  Distance  and  velocity  measurements  for  the  2008  June  1  CME’s  leading  edge 
computed  under  the  “Fixed-0”  approximation,  assuming  two  different  trajectory  angles 
The  implausible  acceleration  far  from  the  Sun  seen  assuming  </>  =  63°  can  be  removed 
by  decreasing  the  trajectory  angle  to  0  =  50°,  though  this  angle  is  inconsistent  with  the 
morphological  analysis. 
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Fig.  4. —  A  demonstration  of  how  the  part  of  a  CME  front  seen  as  the  apparent  leading 
edge  can  change  as  the  CME  expands.  The  dashed  red  line  connects  the  observer  to  the 
true  leading  edge  of  the  schematic  CME  front.  The  green  line  connects  the  observer  to  the 
apparent  leading  edge,  with  the  green  asterisk  indicating  where  the  observer  will  think  the 
leading  edge  is  based  on  the  “Fixed-0”  approximation.  The  apparent  leading  edge  moves 
from  the  far  side  of  the  CME  front  in  the  top  panel  to  the  near  side  in  the  bottom  panel. 
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Fig.  5. —  Our  best-fit  kinematic  model  for  the  leading  edge  of  the  2008  June  1  CME,  based 
on  the  “Harmonic  Mean”  measurements  from  Fig.  2.  The  model  assumes  three  distinct  kine¬ 
matic  phases:  an  initial  phase  of  constant  velocity,  a  second  phase  of  constant  acceleration, 
and  a  final  phase  of  constant  velocity. 


Fig.  6. —  Three  multipanel  figures  comparing  real  and  synthetic  images,  representative  of 
C0R2,  HI1,  and  HI2  data.  Each  frame  consists  of  a  STEREO-A  image  (A),  a  STEREO-B 
image  (D),  our  evolving  3-D  model  CME  (C),  and  synthetic  images  computed  from  this 
model  (B  and  E).  Panel  C  has  a  CME  speed  indicator  in  its  lower  right  corner.  In  the 
HI2  frame,  Panel  C  also  shows  the  locations  of  STEREO-A,  STEREO-B,  Earth,  and  Comet 
Boattini.  Panel  F  shows  the  the  locations  of  the  model  CME,  STEREO-A,  STEREO-B, 
and  Earth  (®)  in  the  ecliptic  plane,  along  with  colored  lines  indicating  the  fields  of  view  of 
the  instruments  used  in  A  and  D.  Panel  G  shows  in  situ  measurements  of  proton  density, 
velocity,  and  magnetic  field  from  STEREO-A  and  STEREO-B.  (A  movie  version  of  this 
figure  is  available  in  the  online  version  of  this  article.) 
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Fig.  7. —  Distance  and  velocity  measurements  for  the  2008  June  1  CME  computed  under 
the  “Fixed-0”  approximation  assuming  a  trajectory  of  <f>  =  63°,  as  in  Fig.  2,  and  also  after 
modifying  the  distances  based  on  CME  front  shapes  suggested  by  our  best  3-D  morphological 
model.  The  modified  velocities  lack  the  implausible  velocity  increase  seen  in  the  unmodified 
measurements. 


